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The detection of the donor star in IYUMa 
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ABSTRACT 

We present the results of a search for the donor star in the high inclination SU UMa 
type cataclysmic variable IY UMa. We detect absorption features in the near infrared 
consistent with an M type dwarf donor star. Using the skew mapping technique to 
exploit the velocity information provided by the 8183-8194A Nal absorption doublet, 
we locate the absorption at the expected donor velocity of IY UMa. 

Key words: stars: novae, cataclysmic variables, stars: individual: IYUMa, stars: 
dwarf novae, techniques: spectroscopic 



1 INTRODUCTION 

IYUMa is a dwarf nova cataclysmic variable (CV). Patter- 
son et al. (2000) (hereafter P2000) and Rolfe, Haswell and 
Patterson (2001) present detailed photometric studies of the 
January 2000 superoutburst and superhump phenomenon. 
Spectroscopic observations have been published, revealing a 
bright hotspot and deep white dwarf absorption lines during 
quiescence (Rolfe, Abbott & Haswell 2001, Patterson et al. 
2000, Rolfe et al. 2002); an eccentric disc during superout- 
burst (Wu et al. 2001); and lines powered by reprocessed 
boundary layer emission during outburst (Rolfe, Abbott & 
Haswell 2002 and Rolfe et al. 2002). 

The orbital period of IYUMa is 1.77 hours, and Patter- 
son et al. (2000) estimated M wd = 0.86±0.11M Q , M donor = 
0.12 ± 0.02Af© and orbital inclination i = 86? 8 ± 1?5. Di- 
rect measurements of orbital parameters are desirable. A 
model-independent solution can be obtained if the orbital 
period and white dwarf eclipse width (already accurately 
measured for IY UMa) can be combined with measurements 
of the donor and white dwarf orbital velocities. The orbital 
period suggests a spectral type of M4.2 to M6.3 for the 
donor in IYUMa (Equation 4 in Smith & Dhillon (1998)). 
Thus absorption features in the far optical and near infrared 
are expected cf. Z Cha and HT Cas (Wade & Home 1988 
and Marsh 1990). We obtained spectroscopic observations 
of IY UMa to look for these absorption features in the wave- 
length range 7000-8300 A. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations were taken on 2001 Jan 3 and 4 using AL- 
FOSC on the Nordic Optical Telescope. A few calibration 
images were taken on Jan 7. A slit width of 1.2" was used 



Table 1. M dwarf template stars 



Object 



Spectral type 



Reference 



GJ 1154, Zkh 176 
GJ 228B, Zkh 84 
BD+28 2110 
CTI 092053.7+280101 



M5V Zakhozhaj (1979) 

M4V Zakhozhaj (1979) 

M3V Upgren (1962) 

M1.5V Kirkpatrick et al. (1994) 



with GRISM 8 yielding a FWHM resolution of ~6.2A and 
dispersion 1.24A per pixel. We discuss only the red end of 
the spectra here (7000-8350A). 39 360-second exposures of 
IYUMa and a nearby comparison were obtained, along with 
exposures of several M type dwarfs for use as template stars 
(listed in Table [l]). Frequent arc lamp exposures provided 
wavelength calibration (giving 9kms _1 r.m.s. velocity vari- 
ation in the sky lines) while HD 93521 served as a flux stan- 
dard and telluric standard. 

Data reduction was carried out using iraf. The CCD 
suffers from fringing but normalized lamp flats were used to 
flat-field the images and removed all detectable traces of the 
fringing pattern. The telluric package in iraf was used to 
correct for telluric features. Apart from the strongest feature 
around 7600A, the telluric correction worked quite well. 

Conditions were not photometric so flux calibration of 
these red spectra was not possible. However, the skew map- 
ping technique relies on normalized spectra. Each spectrum 
was normalized by fitting a 4th order polynomial to the con- 
tinuum regions and dividing by this. 
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Average IY UMa spectra from 3rd & 4th January 2001 
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Figure 1. Average red spectra of IYUMa on 3rd & 4th January 2001. (a) shows average in rest frame of white dwarf, (b) in rest frame 
of donor, (c) in frame of observer, (d) shows average when donor eclipses disc and white dwarf, (e) when disc obscures donor, (f) shows 
M4V template spectrum. Shaded regions show error bars. 



3 AVERAGE SPECTRA 



Various 1 ullTtii'tiii! averaged of Hie IYUMa spectra are shown 
in Figure [i]a.-e along with the M4V dwarf spectrum in Fig- 
ure h|i. Each plotted spectrum has been Gaussian smoothed 
with FWHM of 4 pixels (less than the FWHM spectral res- 
olution), and shifted in flux by a multiple of 0.5 units. 

Figure [l]c shows the average of all 39 IY UMa spectra. 
The most obvious features are the Hei emission and core 
absorption at 7065A, the 7773A O I absorption, the broad 
absorption band around 7600A to 8200A and the absorption 
line around 8200A. The He I is from the disc. 

The broad 7600A-8200A feature matches the TiO ab- 
sorption band seen in the M dwarf spectra, where it is 
strongest in the M5V dwarf. Figures show average 

spectra of IYUMa during eclipse (orbital phases -0.1-0.1) 
and where the donor should be partially obscured by the 
disc (phase 0.42-0.58 assuming a quiescent disc radius of 
0.25a (Stanishev et al. 2001)). The TiO absorption is clearly 
much stronger (by a factor ~2) when the donor eclipses the 
disc than when the disc obscures the donor, confirming the 
identification of this feature as the TiO band absorption in 
the donor. 

The 7773A 1 triplet has been seen in many CVs. A 
study of 65 CVs by Friend et al. (1988) reveals this fea- 
ture in absorption in several systems during outburst, but in 
emission during quiescence. HT Cas also sho ws this feature 
in emission during quiescence ( |Marsh 199c| ). Friend et al. 
(1988) attribute this to line emission from an optically thin 
disc in quiescence, and absorption from an optically thick 



disc in outburs t. Two high inclinat ion systems, Z Cha (Wade 
fc Home 198c]) and V2051 Oph ([Friend et al. 1988;), show 



a broad emission from the Ol triplet, superimposed with a 
narrower absorption component. It is possible that a broad 
weak emission component is also present in these IYUMa 
spectra. Friend et al. (1988) suggest that the narrow absorp- 
tion feature is absorption of white dwarf emission by the 
disc, something which we might expect in high inclination 
systems. That this feature arises from the disc and certainly 
not from the donor is confirmed by Figures [j^, & b, where 
each individual spectrum was Doppler-shifted into the rest 
frame of the white dwarf or donor before adding together 
the spectra. The velocities of the donor and white dwarf 
were assumed to be Ka = 442 kms _1 and K\ — 62 kms -1 
respectively, calculated using the P2000 orbital parameters. 
In the rest frame of the white dwarf the O I feature remains 
sharp, while in the donor star frame it splits into two ab- 
sorption dips. In addition, the depth of this absorption is 
significantly stronger when the disc is in front of the donor 
(Figure ^e) than when the donor eclipses the disc (Figure 

W)- 

The 8200A absorption feature closely matches the 8183- 
8194A Nal doublet seen in all four dwarf template spectra. 
In the average spectrum (Figure ^c) it does not have the 
same doublet structure as seen in the template spectra, but 
when shifted to the donor rest frame (Figure |l|b) the doublet 
structure becomes apparent, clearly matching the laboratory 
wavelengths. This feature is also more significant (consider- 
ing the flux errors) during disc eclipse (Figure [l]d) than when 
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the disc passes in front of the donor (Figure |l|e). This is a 
clear detection of the Nal doublet, a feature which has been 
exploited to measure Ki in other CVs, c.f. Wade & Home 
(1988), Friend et al. (1988, 1990a, 1990b) and Marsh (1990). 

The TiO absorption band around 7000A-7400A is also 
more prominent when the donor obscures the disc than when 
the disc obscures the donor. The K I lines at 7665A and 
7699A are not clearly seen, but there are weak features close 
to the noise level which could be tentatively identified as 
these lines. 

Weaker Nal and TiO absorption around phase 0.5 than 
at other phases has been see n before in dwarf novae , e.g. 
Z Cha ( |Wade fc Home 198Sj ), HT Cas ( |Marsh 199(j ) and 
SS Cyg ( |He^srn^n*^t*lil**T984 ) , and was attributed to a non- 
uniform distribution of absorption across the surface of the 
donor. The suggested cause of this non-uniform distribution 
is the filling-in of absorption features on the inner face of 
the donor by irradiation from the disc and/or white dwarf. 
This, rather than eclipse of the donor by the disc, might also 
be the cause of the weakened features around phase 0.5 in 
IYUMa. 

The absorption feature at 7600A is probably a com- 
bination of residual atmospheric absorption and the 7600A 
feature seen in M dwarf spectra. There is some evidence of 
over-correction of the telluric absorption in the M5V tem- 
plate around 8230 A. 

We see no evidence of Paschen emission around 8200A 



seen in some CVs e.g. TT Ari and V603 Aql (Friend et al 
198SjT| 

The donor star in IYUMa has been detected, showing 
all the expected features of a late M dwarf star. Without flux 
calibrated spectra, the relative strengths of the absorption 
features cannot be used to estimate more precisely the spec- 
tral type of the donor. The low signal-to-noise ratio prevents 
us from exploiting the Nal doublet to measure Ki directly 
f rom the spectra, so we adopt the s kew mapping technique 
(smith, Cameron & Tucknott 1993) to locate the donor in 



velocity space. 



4 SKEW MAPPING 

Skew mapping is the best method for detecting the donor 
and determining its velocity when the absorption features 
are weak. This technique has been successfully employed 
for several systems (e.g. Smith, Cameron and Tucknott 
(1993), Smith, Dhillon and Marsh (1998) and recently in VY 
Aqr (Littlefair et al 2000)). Skew mapping co-adds spectra 
to increase signal-to-noise, while also correcting each spec- 
trum for the orbital motion of the donor, to avoid smearing 
out donor absorption features. Skew mapping extends the 
technique used in producing Figure [ljh> (where K% was as- 
sumed), by cross-correlating each such co-added spectrum 
with a template donor star spectrum, and finding the com- 
bination of template star and donor velocity at which the 
cross-correlation peak is strongest. This process is equiv- 
alent to producing a trailed "cross-correlation spectrum", 
and then back-projecting this trail into velocity space, pro- 
ducing a cross-correlation map in velocity-space (the skew 
map) which should show a peak at the velocity of the donor 
star. 

The normalized spectra of the template M dwarfs and 



IYUMa were continuum subtracted, and were Gaussian 
smoothed with FWHM of 2.5 pixels to increase S/N with- 
out losing any spectral information. Each IYUMa spec- 
trum, F(X), was cross correlated with the template spec- 
trum, T(X), for a range of velocity shifts, v, using correlation 
coefficient C v defined as 



C v — 



£A F(A)T ' u(A) 



(i) 



where T' V (X) = T(^j^). This definition of T' V (X) simply 
has the effect of red-shifting the template spectrum by ve- 
locity v before cross-correlating. The denominator of C v is 
chosen so that if the template and donor spectra are identi- 
cal, the value of C v at the peak will be 1, independent of any 
simple scaling of the flux of either spectrum. The correlation 
trails thus produced were then transformed to velocity-space 
using the Fourier-filtered back-projection technique as im- 
plemented in Tom Marsh's molly software. The strength of 
the skew map at a given velocity therefore shows how well 
the object spectra match a star of the template's spectral 
type orbiting at that velocity. 

Only the Nal doublet wavelength range was used in 
the cross-correlation, since this is the only sharp donor ab- 
sorption feature detected, and will therefore provide all the 
velocity information. Skew maps were also produced us- 
ing the spectral range A = 7045-8240A (with the O I and 
He 1 lines masked out) which uses all the detected donor 
star features, but this simply smears out the donor star in 
the skew map, probably a result of the broad TiO absorp- 
tions showing a strong correlation but low velocity sensitiv- 
ity. The normalization of the spectrum around the Na I dou- 
blet was improved by fitting a straight line to the contin- 
uum either side of the doublet. The skew map also depends 
upon the relative systemic velocity of IY UMa and the tem- 
plate star. Neither of these is yet known, although measure- 
ments of Ha emission lead to estimates for IY UMa of 7 = 
lS.Oilkms-^Rolfe et al. 2002) and 7 = -4±32 kms" 1 (Wu 
et al. 2001). Such estimates should be treated with caution, 
however, since the determination of systemic velocity from 
emission lines is notoriously unreliable: effects such as phase- 
dependent occultation/absorption of emission and tempo- 
ral variations in emission can produce incorrect results. By 
cross-correlatin g the M3V t emplate whose systemic veloc- 
ity is lkms" 1 ( Evans 1967 ) with the other four templates, 
the systemic velocities of these stars were determined. The 
skew maps for each template were produced for a range of 
values of the IYUMa systemic velocity, 7iYUMa, between - 
100 and 300 km s -1 (the low resolution of our spectra limits 
the accuracy of our systemic velocity determinations). This 
was identified as a suitable range for 7iYUMa after initially 
studying the range -750-750kms _1 . This procedure enables 
a determination of the value which produces the clearest 
peak. We do not consider this a determination of the sys- 
temic velocity of IY UMa. 

The orbital phase range 0.42 to 0.58 was not included 
when producing the back projections to avoid the weak- 
ened/obscured absorption lines around this phase interfer- 
ing with the results. The results from using the entire phase 
range are considered at the end of this section. 

Figure |^ shows a phase- folded and -binned and velocity- 
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Figure 2. Cross-correlation trail for the M4V template 
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Figure 3. IY UMa skew maps using the M4V template. The 
teardrop shows the expected velocity of the donor, and the black 
crosses show the white dwarf, system and donor centre of mass 
velocities. The white cross is the measured peak maximum. 



binned correlation trail for the M4V template. All four tem- 
plates yield almost identical trails; at this low resolution and 
signal-to-noise, differences in structure in the Nal doublet 
between the spectral types are not going to be distinguish- 
able in the IYUMa observations. The trail shows an S-wave, 
corresponding to a localized region in velocity space at which 
the IY UMa spectrum best matches the sodium feature. The 
peak-peak amplitude of this feature is about 1000 km s , 
and is at minimum excursion around orbital phase 0. This 
phasing and amplitude is in agreement with the predicted 
orbital motion of the donor (plotted in white in Figure ^ 
using the P2000 orbital parameters (line-of-sight donor ve- 
locity semi-amplitude 442 kms -1 ). The velocities are helio- 
centric and corrected for the template systemic velocity. 

The trails were transformed to velocity space. The re- 
sulting skew maps have a clear maximum close to the pre- 
dicted location of the donor star using the P2000 param- 
eters. The exact location of the maximum depends on the 
assumed IYUMa systemic velocity, 7iYUMa- Two skew maps 
for the M4V template are shown in Figure []. The map in 
the left panel uses the value of 7iYUMa which produces the 
strongest peak (i.e. max correlation), while that on the right 
corresponds to the 7iYUMa which produces a peak with no 
x-component of velocity (as expected for the donor) . These 
values of 7iYUMa were determined as described below. 

The velocity and width of the peak was determined for 
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Figure 4. Comparison of skew maps for each template. 



each image by fitting a 2-dimensional Gaussian of the form 



CV = Co + Ci exp 



u x .max 



2Av 2 r 



2Av% 



where CV is the value of the skew map at velocity (v x ,v y ). 
Figure ^] shows the variation in position and shape of the 
skew map peak as a function of 7iYUMa- The phase shift 
is the angular offset of the peak from v x =0 measured about 

— . The velocity ampli- 

The peak width is measured 



the centre of mass, i.e. ^- arctan 

7 Zty 



tude is V = yj vl, max + 



using Av r = sj Avx + Avy. The correlation of the peak is 
taken as Co + Ci . 

As 7iYUMa increases from -100 kms -1 to 300 kms -1 , the 
peak moves from a velocity of about 600kms -1 and phase 
shift of around -0.01 through the expected donor location 
to a lower velocity around 400-450 kms -1 and phase shift of 
around 0.13. This result is true for for all four templates. 
Only the peak size A?v and the peak correlation vary ap- 
preciably depending on the template used. The difference 
in Av r between templates should not be taken as a mea- 
sure of how tightly each template constrains the location of 
the donor: it is more likely to result from differing seeing 
between the different template observations. However, the 
variation of Av r with 7iYUMa for each template is more use- 
ful, telling us which value of the IYUMa systemic velocity 
gives us the sharpest peak, and is therefore the most likely 
value for the systemic velocity. This shows the peak to be 
narrowest for 7iYUMa around 75 to 125 km s -1 . The most im- 
portant measure of which value of 7iYUMa is correct is that 
for which the peak correlation is strongest. This is around 
7iYUMa=100 to 150 kms -1 for the four templates. The peak 
correlation of the templates is similar, which considering the 
low signal-noise and resolution of these observations is not 
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surprising: detailed differences in the Nal doublet between 
the spectral types have not been resolved. 

The value of TiYUMa with the highest peak correlation 
gives us one measure of the best fitting skew map. The cor- 
responding M4V skew map is shown in the left panel of 
Figure |§| The average phase shift for all four of these maps 
is 0.044, with a spread in values of 0.001. The 360 second ex- 
posure time used corresponds to a phase resolution of 0.06, 
so this phase shift is consistent with the peak actually lying 
along the line of centres of the two stars, where we expect 
the donor. The accumulated error in the P2000 ephemeris at 
the time of these observations is only 0.005 orbits. This then 
gives us a donor velocity of iG=383 km s _1 with a spread 
of values of 6kms _ . If this measured phase shift is real, 
then we should also consider values of 7iYUMa for which the 
phase shift is zero, placing the donor along the line of cen- 
tres. The corresponding skew maps (shown for M4V in the 
right panel of Figure pi) give iG=524 km s _1 with a spread 
of 19kms -1 . A possible cause of such a phase shift would 
be irradiation of the donor by emission from the stream-disc 
impact, filling- in the Na absorption mainly on the leading 
side of the donor, thus moving the correlation peak in the 
direction seen. 

As described earlier, these skew maps were produced 
omitting phase range 0.42-0.58. As a check, the analysis 
was repeated using the entire phase range. The skew maps 
corresponding to phase shift zero show far less distinct peaks 
than in Figure ^, rendering this measurement less reliable. 
The peak in correlation, the minimum value of Av r and 
phase shift of zero occur for values of 7iYUMa about 10- 
20 kms -1 greater than before. The phase shift at maximum 
correlation is now 0.057 with a spread of 0.003, as large as 
the phase resolution of the data. The orbital velocity of the 
donor for maximum correlation is now K2—362 kms -1 with 
a spread of lkms -1 . 



5 CONCLUSIONS 

The donor star has been unambiguously detected in IY UMa. 
It is a late type M dwarf, but the low signal-to-noise and 
resolution of these observations prevents an accurate deter- 
mination of the spectral type. The donor velocity, K2, lies in 
the range about 380-540 kms -1 , consistent with the value of 
442 km s -1 using the model-dependent orbital parameters 
estimated from the photometric study in P2000. Observa- 
tions with a larger telescope, providing higher spectral and 
temporal resolution and signal-to-noise would facilitate a di- 
rect and accurate determination of both the spectral type 
and velocity of the donor star in IY UMa. 
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